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Abstract 

An exact macroscopic extended model for ultrarelativistic gases, 
with an arbitrary number of moments, is present in the literature. 
Here we exploit equations determining wave speeds for that model. 
We find interesting results; for example, the whole system for their 
determination can be divided into independent subsystems and some, 
but not all, wave speeds are expressed by rational numbers. Moreover, 
the extraordinary property that these wave speeds for the macroscopic 
model are the same of those in the kinetic model, is proved. 



1 Introduction 

The macroscopic extended model, with an arbitrary number of moments, for 
ultrarelativistic gases has been introduced in [1] which is the generalization 
of [2]. It proposes the field equations 

= r^-"", for n = 0,l,...,A^, (1.1) 

where A"'"^"'"'" and are symmetric and trace-less tensors. In [1] it is 

proved that the entropy principle for this system amounts to assuming the 
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existence of the symmetric and trace-less Lagrange multipliers Xaia2---a„ and 
of an arbitrary function F{xq ), such that 

A"""-"™ = j F^{X, AX, ■ ■ ■ , K.-mP^' ■ ■ -p^npy^' ■ ■■p'^'-dP, (1.2) 

1 n 7 7-. dp^dp^dp^ , , . 

where rm = ttt^, dP = and is the four momentum satis- 

fying the relation p^p^ = (because the particle mass is zero in an ultrarel- 
ativistic gas). Moreover, Xai---a„ can be taken as independent variables. A 
further condition has been studied in [3] , while a reader not acquainted with 
the general context of Extended Thermodynamics might profit by reading 

mm. 

In this article we study the wave speeds for the system (11. ip . To this end 
we observe that this system, taking into account (II. 2p . can be written as 



AT 



J] A""i-"'"'^^-^"a„A/3,...;3„ = ri-"'^, (1.3) 



n=0 

with 



=j F,n,„(A, A^, ■ ■ ■ , A^,...^^p^i ■ ■ ■p^'-)pVi ■ ..p'^-pP- . ..pf-dP, (1.4) 
where Fm,n = ^ „ . Now recall that hyperbolicity for the system (11.31) 



in the time-like direction (with ^q,^° = —1), in the sense of Friedrichs 
[3 [HE], means that 

1. The system ^a^""^'""'"^^'"'^"5A^i...^„ = 0, in the unknown (iA_a^...^„, has 
only the solution dXp^.-.p^ = 0; 

2. The system 

^^A^-i-^-^i-P"5X^,...ls^ = 0, (1.5) 

where (pa = Va~ X^a, has real eigenvalues A and a basis of eigenvectors 
'^A/3j.../3„, for every unit vector r]a such that = 1, fjaS." = 0. 
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The eigenvalues A are the wave speeds to obtain. A sufficient condition for 
hyperbohcity in every time-hke direction can be found in [TU]. The hyper- 
bohcity for our system (11 .Sp has aheady been proved in [Tj as a consequence 
of the so-called "convexity requirement" ; so it remains to find the values of 
the eigenvalues. We ffist introduce the time- like direction Ua = 7~^Aa with 
7 = y/— A^A^ and the projector 

Kp = gap + UaUfi (1.6) 

into the 3-dimensional subspace orthogonal to u^- After that we note that 

h"-^ipi3 7^ 0, otherwise = —(p^u^Ua and the condition 2. becomes 1. with 
= Ua, from which dA^^...^.^ = because of the fact that d\p^...i3^ is an 

eigenvector. Because h°'l^ipp ^ 0, we can now define Lp = ^JhP^^hp^^ and 
= V'^^aVp, so that Va is a unit vector orthogonal to Ua- Finally, we 

introduce Kafi = ha/s — VaVjs which is the projector onto the 2-dimensional 

subspace orthogonal to Ua and Va- 
in [H] we have found some properties satisfied by these projectors and 

the definition of K^^ ■ ■ ■ K'^'^ , a tensor such that for every tensor T^^"'^p 

the tensor T<^ ^i->> = ■ ■ ■ Kf^^Tl^^-l^^ is equal to the sum of the 

expression Kj^^ ■ ■ ■ Kj^T^'^"''^p and a linear combination of its 2-dimensional 
traces (i.e., i^^,^, ■ ■ ■ ■ ■ i^jT^i---^-) and, moreover, the 2- 

dimensional trace of T^^ 7i --7p> jg equal to zero. In other words, T^^ 7i- -7p> 
is the 2-dimensional trace-less part of Kj^^ ■ ■ ■ Kj^T^^"'^p . In Section|2]we will 
report some of these results. We recall that, in [11], only the case = 3 has 
been treated, while in the present paper we study the general case. However, 
some details of [IT] are also general, so we report these results here, for the 
sake of completeness. Obviously, the objective of all these considerations is 
to rewrite the system (II. 5p in terms of tensors belonging to our 2-dimensional 
subspace. This result will be given in Section |2l 

The system (II. 5p will not assume a very elegant form, a consequence 
of its generality. For this reason in section |3] we will consider the cases 
p = N,p = N— 1 and p = N — 2. In this way the case N = 2 will be 
exhausted. The resulting wave speed, in the reference frame moving along 
with the fluid, are A = 0,A = |, A = ±l, A = ±^. These are expressed by 
rational numbers, except for the last one. 

Finally, in Section H] we will prove the extraordinary property that the 
wave speeds do not depend on the function F, so that we can take F = 
Q-T,i^o^^/^ where K is the Boltzmann constant; in other words, the wave 
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speeds for the macroscopic model are the same as those in the kinetic model 



2 On the 2-dimensional trace-less part of a 
tensor and other properties of the projec- 
tors 

Let us begin by exploiting some technical tensorial properties which will be 
useful to write the system (11 .Sp in terms of 2-dimensional traceless tensors. 
First, let us define the tensor 

[f] 

= E «^ ^(^^f^^ ■ ■ ■ Kp...,p.K^s::: ■ ■ ■ K)^'''' ■ ■ ■ ^'""^'"^ ' 



s=0 



(2.1) 



where = {—j)^ {p-2s)\ ^^"^r^'" ^^'^ from now on the square brackets denote 
the integer parts of the number. 

What does this tensor mean? When we contract it with a generic symmet- 
ric tensor T^^-Pv we obtain that T<^ ti->> = K^^ ■ ■ ■ K>>T'^i-/3p equals 

the sum of T^^"'^pKJ^^ ■ ■ ■ Kj'' (because oq = 1) and of a linear combination of 
the 2-dimensional traces of T"'^'""'^ . Moreover, the following theorem holds: 

Theorem 2.1 The 2-dimensional trace of T^^ "'^'""'^ is zero, or, equivalently 

K<y^ . . . K,,,, = . (2.2) 



This theorem has been proved in |TT] the other theorems of this sec- 

tion. By means of Theorem 12.11 it is now natural to call T^^ 7i- ■■7p> ^j^e 
2-dimensional trace-less part of T''^i ' >. Then (12. ip gives the 2-dimensional 
trace-less part of xj^^^ ■ ■ ■ Kj^^; it will be useful for the sequel to have a sort 
of inversion of this relation. It is provided by 

Theorem 2.2 The following identity holds: 

[i] 

s=0 

(2.3) 
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with 



r! 1 (2r-4.)!! ^^.4) 



' (r-s)! (2s)!! (2r-2s)!! 

Also the proof of this theorem is reported in pjj. We conclude this list of 
theorems with 

Theorem 2.3 The following identity holds: 

T\r<^ 71 . . . /^7p> /^{"1Q!2 . . . 7>'02s-l«2sT/a2s+l . . . T^a2s+c r ra2s+c+i . . . Tra2s+c+d) 
ai Op 

^ (2g)!! (^■^ + + ^ ~ P)' jy<2 71 T<r^P> ^f°l °P+l /<f"2 «p+2 T^apa2p 

(2s-2p)!! {2s + c + d)\ 

_ J^ Ct2p+ia2p+2 _ _ _ J^ a2s-ia2s '^/ 0'2s + l . . . '^/ a2s+c Jj a2s+c+l . . . JJ'<^2s+c+d_ ^2 5^ 

where the indices underlined denote a symmetrization over those indices. See 
[TT] for the proof of this theorem. 

In [11] we have also proved that the equation 

= ^Z"U,p,. ■ ■ ■ U,p,,V,^,,,, . . . y,,,,,,dA^--^" , (2.6) 

represents an invertible transformation from the variables d\p^...p^ (with zero 
4-dimensional trace) to X'J^^ ^j^^ (belonging to the subspace orthogonal to f/y^ 
and V^^ and with zero 2-dimensional trace) and the condition that dXjs-^.../^^ 
has zero 4-dimensional trace is "translated" into the fact that ^^J^ has 2- 
dimensional zero trace. We note that the index h denotes how many contrac- 
tions with U... are involved, while the index k denotes how many contractions 
with V... are present. 

Therefore, we can now consider X'^^^ ^J*" as unknowns, instead of d\p-^...i3^. 
Obviously, the invertibility above mentioned holds for any symmetric ten- 
sor dXai-.-an with zero 4-dimensional trace; consequently, it holds also with 
LpaA°'^^"'^"^°'^"'°'"- dXa^.-.a^^ iustcad of dXai---a„- othcr words, we have that if 

N 
n=0 



(2.7) 



for all a, b such that a + b < m, then also (II. 5p will be satisfied. The converse 
is trivial. Consequently, (12. 7p is the new system to be used to obtain the 
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unknowns Xj^'J^^' But it is better to rewrite (12 .7^ as 



N 

E i-^.U^U^ + ^Vo) K<^ Kl >Up^_^, ■ ■ ■ U,^^V^^^^^, ■ ■ ■ V,^^^^, 



n=0 



■ v4°''i-^™"i-""dA«,...„„ = , (2.8) 

with p = 0, . . . , m; a + b = m — p. 

It is interesting to see that, from [TT], it is not necessary to substitute the 
relation hnking dXai-a„ to X'^^ J^"^ in this equation, since it will be a natural 
consequence of the kinetic expression (ll.4p . The effective steps can be found 
in [TT], so we limit ourselves to report the result here, i.e., 

[m + n — a] r2r — bl 

io ^-'"^^'^^^^ ^{2r-2s-by. {m + n- ay. (2s - 2p)!! 

■ 1 rp j y ^) ^(^n-p-2r+b+2s-2T,2r-b-2s+2Ty 

N [ 2 J [ 2 J S-p / , , 1 \ 

F (X n.)n-^-^'dn (^r - 2.)! (m + n-g-fe-p)! (2.)!! 
Pm,n[^,Pl)P ""P ^2r-2s-b-l)\ (m + n+l-ay. (2s - 2p)!! 

2^ j V^^; "''■(n-p-2'r+6+2s-2T+l,2r-6-2s+2T-l) 

To state the system (12. 9p we have to choose a value of p such that < 
p < N; after that we have to write (12. 9p for every m such that 

p<m<N (2.10) 

and for every couple (a, b) such that a + 6 = m — p. In this way a set of 
equations is obtained which are linear combinations of tensors X'^^j^^'', all 
of the same order p; moreover, the coefficients of these linear combinations 
are scalar functions. In other words, the wave speeds are obtained by simply 



m+71+l— a 



V ^/ ^(n--D-2r+b+2s~2T+1.2r-b-2s+2T-l) ^ ' K"^-^) 
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imposing that the matrix of these scalar coefficients is singular. Doing this 
for every value of p, we obtain all the wave speeds. As a result we have 
proved that the whole system, for the determination of the wave speeds, can 
be divided into some independent subsystems for each given value of p. We 
note also that for every fixed value of p, from (11. 5p only those with m > p 
have to be considered in the system (12.91) : moreover, in this system we have 
n > p so that, from the initial unknowns ^A^^.../^^ only those with n > p are 
present in the system (12. 9p . 

3 The wave speeds for some values of p 

3.1 The wave speeds for p = N 

In this case, from (I2.10p . we have m = A^. In other words, only the last 
equation of system (11.51) contributes to the system (12.90 . Moreover, a + b = 
m — p imphes a = 6 = 0. As a consequence, we see that (12.91) becomes 

N 

because, in the coefficient of ip the summation appears, so that this co- 

n=N+l 

efficient is zero. Moreover, we have defined Gm,n = /o°° -^m,n(A, pj)p'^^"''^^dp. 
Consequently, we have (p^W^ = as unique eigenvalue. It corresponds to the 
waves moving along with the fluid, because in the case = [/^ we have 
= (fi^U^ = {—XU^ + ri^)U^ = A. In other words, the wave speed is zero. 

We also note that, ip^U^ = Q does not depend on Fm,n, so it is the same 
as for the kinetic model |12] . 

3.2 The wave speeds for p = N — 1 

From (I2.10p we obtain m = A^ — 1 or m = A^. In other words, only the last 
two equations of the system (II. 5p contribute to the system (12. 9p . 

If m = A^ — 1, from a + h = m — pii follows a = 0, 6 = 0. If m = A^, from 
the same equation we obtain a + 6 = 1, i.e., a = 1, 6 = or a = 0, 6 = 1. In 
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this way we obtain three equations, 
4-7r 

~ '^^^ \2N - ~ ^^-^^^ ' ^^-^'^-1 ^(0,0)'"" + ^^-1.^ (3-1) 

~ ^/-^' (2iV_l)! (^ - l)!(2iV - 2)!! Xj^^^- + G^,^ (3.2) 

■ [^Mf/'^^^]^iV!(2iV - 2)!!X--- + ^jJ^NNm - 2)!!X---] = , 



. r _ 87riV(iV)! _ 2)!!j^7i^-7iv-i _ jy/' ^, iViV!(2iV - 2)!!XZ;-,-;^^-^l = 

L ^(2iV+l)!^ ' ^'^ (2X+1)! ^ ' ^ 

respectively. Now we note that (13. ip and (13. 2p can be considered as two 
equations in the two unknowns 

- ^^f^" (^Ar^(^ - l)'(2A^ - 2)!! and 

^.f^^ (^^^Ar!(2iV - 2)!! Xj-^- + V^^^^^iViV!(2iV - 2)!! Xj--- , 
where 

is the coefficient matrix. But, in the particular case v? = 0, we would obtain 
only the zero solution for the equations (I3.ip - (l3.2p and this for condition 1. 
of hyperbolicity; therefore the matrix (13. 4p is nonsingular and Gjy^N 0, so 
that the two unknowns are zero, that is 



N 



y,,U" Xjif- + 2^ _^__^A-;-™-' = . (3.5b) 
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Also fl3.3p has to be imposed. Now, if ^^^U^ = 0, these three equations 
become 

(3.6a) 

G^.^X^-^;--'^-- =0, (3.6b) 

so that there is a free unknown. Consequently, ^p^U^ = is an eigenvalue. 
If 7^ 0, (13. 5p and (13. 3p can be written as 

j^7lv,7iV-l _ g 

so that the eigenvalues have to be determined from the equation 

Thus (13. 7p does not depend on Gm,n and hence it is the same as in the 
kinetic case. Instead, the eigenvectors defined by (13. 6p depend on Gm,n, so 
that they are not the same as in the kinetic case. In particular, if = ?7^, 
from which ip^U^ = \ and Lp = 1, the wave speeds are A = i^^^^ which are 
rational numbers satisfying |A| < 1, which means that the wave velocity does 
not exceed the speed of light, as expected. Moreover, our results is consistent 
with [13]. 



G 



71. 
(0,1) 



■,7iv-i 



0, 



Y'- 



7i,...,7]V-i 



3.3 The wave speeds for p = N — 2 

This part is important because, if = 2, it completes the set of wave speeds 
for the 14 moments model. From (12.101) we now have m = N — 2 or m = N — 1 
or m = N] more precisely, we have to write the (12. 9p with 



• 


m = 


N - 


2, 


a = 


0, 


b = 





• 


m = 


N- 


1, 


a = 


1, 


b = 





• 


m = 


N, 




a = 


2, 


b = 





• 


m = 


N- 


1, 


a = 


0, 


b = 


1. 
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• m — N, a — 1, b — 1, 

m m = N, a = 0, 6 = 2. 

In this way we obtain the system 

-iA7i---7iV-2 I ^ ■iA7l---7iV-2 



2 +Gn-2,nX3 ^— 

x^71---7iV-2 , ^ ■v^7l---7iV-2 I ■v^7l---7iV-2 n 

^N-l,N-2-^l +^N-1,N-1^2 +(j7V-l,Af^3 — U, 

vTi---7iV-2 I r-i vTi---T^-2 _ n 



'7i---7iV-2 _|_ ^ 

^N,N-2^l + 

-1A71---7JV-2 , ^ 
LrAr-l,iV-2 J^l + 

TA7i---7Af-2 I V' 



\rll---lN-2 , -1A71---7JV-2 _ ri 

Ar-l,JV-l-i2 + tJiV-l,iV -1^3 — U. 



7l---7Af-2 , -1^71 ■■ -7^-2 n 

+ ^N, N I3 — U , 

Gryl^—lN-l I ry-^-L—lN-l , rv^\—lN-1 rv 



where 



An 



{N - 2)\{2N - A)U Xj^--^:^''-' 



{2N - 3) 
An 

^7i...7iv-2 ^ ^^[/M__^!L_(Ar - l)!(2A^ - 4)!!X,^^-^^-^ 



(2A^-3)! 



An 



+ ^ (2iV-l)! ^^^ " ^^^^ " ^^-^^^ " 



(1,0) 
71---7JV-2 



j^7i---7jv-2 



An 1 

(2Ar-3)!2 ^ ' (2,0) 



4)!! x?'"'^'^"^ 




{2N-1)\2 
An 1 



N\{2N - 2)11 Xj^--^:^"-' 




, -X2N -\){N - 1)2- ^(27V - 4)!! X^,}":. 

(2iV + l)!^ ' (2A^-1)!^ ' (I'l) 



71---7JV-2 



71---7JV-2 



An 



7i...7iV-2 



- <P 



^(2F3l)!2'V(W-l)!PJV-4)!!X- 

An 



(3.8a) 
(3.8b) 
(3.8c) 
(3.8d) 
(3.8e) 
(3.8f) 



(3.9a) 



(3.9b) 



(27V- 1)! 



(]V-l)2(iV-l)!(2iV-4)!!X(J;) 



7i---7iV-2 



An 



{2N 



2{N -1){N- 1)\{2N - 4)!! Xj^'-^"-' 



(3.9c) 
(3.9d) 

(3.9e) 



10 



+ 6(/7 ■ 



4:71 N\ 

Air 
An 



''(2]V3T)!('^-1)^K2A'-4)!!A7';„; 



JV(JV - l)]V!(2iV - 4)!! Xj' „ 



71...7JV-2 



47r Att 



(2Ar + 

+ (^]^^^'(2iV - 2)!! X--- - (^^A^A^!(2iV - 2)!! X^' 

(3.9f) 

Z?-'"-' = -^,U^ (^^^!^), 2(iV - l)!(2iV - 4)!! Xj-^-- , (3.9g) 
Z^^---- = -V^.^" (^^^^^), 2(iV - l)(iV - l)!(2iV - 4)!! X^-^-^^- 



+ ^ pr^l2(iV - l)iV!(2iV - 4)!! Xj-^-- , (3.9h) 

A-K 



~ ^^^' (2iV^Il)! ^^'^^^ " ^^^^^ ~ ^("' 
47r 

^ --7i...7jv_2 



+ V^M^^^ (2^^^^' ,„ iViV!(2iV - 2)!! X^--^-- 
47r 



^ (2iV + l)! ^^^^^ ~ ^^^'^^^ " ^S-?"" • ^^-^'^ 

Now, (l3SiD-(l33c]) give x7^-^^-^ = 0, X^'-"""-' = and Xg^^-^^"^ = 0. 

If = 0, fl3.9ap becomes an identity so that we have five equations in 

the six unknowns Xj;,-;--, Xj-'^--, X^-^^^;--, Xj-^-/--, X^-^-'/--, and 
"^(0 2/^"^' consequently ip^lf^ = gives a wave speed. 

if yj^f/'^ ^ 0, from we have XJ q)'^'^-' = 0. From (IHldl), (^^9^, 

and dHh]) we obtain y^^-^^-^ = 0, and Z^'-^'^-^ = O; after that (ISIM 
and fl3:8e|) give y^^^-^^-^ ^ g ^^^^ y^7i-7iv-2 ^ Consequently, we have the 
equations 

^7l...7iV-2 ^ g ^ J^7l...7iV-2 ^ g ^ y7l...7iV-2 ^ g ^ y7l...7iV-2 ^ g ^ 
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together with to determine Xj^-^""-', X^^l-^""-', 2)''"'-', X^to^''^ 

^Jir"'- Now X^'-^''-' = and Y^'-^''-' = is a subsystem of two 
equations in the unknowns Xj^''^^''-^ X^^^f''-\ If ((/J^t/'^)^ - ^ip^ = 0, 
one of these equations is a consequence of the other, so that we have four 
equations to determine five unknowns. Therefore, we have obtained another 
wave speed. Also in this case it is smaller than the speed of light, but it is 
not a rational number. 

If {ip^U^)"^ — 2N^i 'P^ 7^ 0' then the equations 

^7i...7iv~2 _ y^"'''""^'' = 0, 

give Xll;^^''-' = and Xj^)''"'-' = 0. From (EUD it follows that zl'-^''-' = 
0, so that (13.8f[) becomes z'l^'"'^'^''^ = 0. Now we are left with the equations 

to determine the unknowns X'^^'^J^''^ , X^'^^^j^^"^, X'^^"^'^^'^ . So the last wave 
speeds are obtained when the coefficent matrix is singular, i.e., when 

0, 



-2N{2N + l)ipf,U'' -{2N-l)ip 
if {2N + 2)ip {2N-l)ip^U^' 

-if^U^' -{2N + 2)ip^Uf' -Sip 



where we have dropped some factors. The solutions are ip^W^ = and 

i^.un' - = 0. 

To conclude this subsection we report the wave speeds found, i.e., 

^,u^ = o, i^,un' = ^^^', = ^uh~y ' 

where 2N-1 — 2N-1 — ^ ^^^^ ^^^^ '^^^ inequality holds only for > 2. On 
the the other hand, if < 2, the case p = N — 2 has not to be considered). 

In the particular case N = 2, which corresponds to the 14 momements 
model, we have found all the wave speeds, namely 

^,u^ = o, ip>,un' = ^v^ {ip,un' = l^\ i^,un' = ^^ 

where the last equation gives the speed of light. 

We note that, in any case, we have found something like 

{'P^U^'Y = kp>^, with < A; < 1. (3.10) 
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On the other hand, the wave speed A is defined by 

[(-A^/. + V,)UnT = kh^^i-K, + V,){-Kv + V.) , (3.11) 

as in the definition of hyperbohcity reported above. Obviously, if = U^, 
then (13. lip becomes = k; in other words, ±\/k is the wave speed in the 
reference frame moving along with the fiuid. So, the following question arises 
immediately: If < /c < 1, will A satisfy also the condition |A| < 1? The 
answer is affirmative and this will be proved in the next subsection. 

3.4 The wave speed in every time-like direction 

The aim of this subsection is to prove that, for < A; < 1, the solutions A of 
( 13. lip will be such that — 1 < A < 1 for every 77^ satisfying the conditions 

^^^ = -1, ^^^'^ = 0, r/X = l. 

Let us consider the reference frame where = (1, 0, 0, 0) and 77^ = (0, 1, 0, 0). 
In this frame fl3.1ip can be written as 

(-Af/° + U')^ = k{l - A') + A:(-Af/° + U')\ 

i.e., the wave speeds A are the solutions of /(A) = 0, with 

/(A) = A' [{Uy{l -k) + k]- 2\U''U\l -k) + {U^f{l -k)-k, 

where we note that the coefficient of A^ is positive because < A; < 1. The 
first question is: Are the roots of /(A) real? The answer is affirmative because 

j = k{l- k) [{U^f - [U^f] + e = k + k{l-k) [([/2)2 ^ (^3)2j > 

3 

where we have used the property — ([/°)^ + ^^(^7*)^ = —1. Moreover, we 

i=l 

easily see that 

f{±l) = il-k){U'TU'f>0, 

i/'(±l) = ±i(f/°)^(l -k)± ^-{Uni -k)- U'U\l -k)±k. 
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Subsituting {U^f = 1 + E?=i(^')^ in the second term with {U^y, it follows 
that 

lfi±l) = ±^(1 - k){U' T U'f ± \ [{U'f + [U'f] (1 - fc) ± i(l + fc), 

so that /'(I) > and /'(-I) < 0. From /(I) > and /'(I) > we obtain 
A < 1, while from /(—I) > 0, and /'(—I) < we obtain A > —1, as expected. 



4 Independence of the wave speeds on ^^ 

Let us now prove that the wave speeds do not depend on Gm,n and, hence, 
on Fm,n- This fact implies that they are the same as for the kinetic model, 
where 

1 poo -| 

Fm,n = ^e-^(^+^-^'^), and = j _e-i(^+^^)p™+"+2^p. (4.1a) 

To obtain this result, let us rewrite (12. 9 p with a = m ~ p — b. We note that 
now the index m is present only in Gm,n- Therefore, (12. 9p can be written as 

N 

J2Gn.,nYJX-''' = 0, (4.2) 

n=p 

for b = 0, - ■ ■ , N — p, m = p + b, ■ ■ ■ ,N. Here, by the definition of the tensor 
'Ybn ''^''y some special cases have to be distinguished. In particular, 

• case 1: b is even and n = p. YfJl^'"'^'' is defined by 
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case 2: 6 is even and n — p + 1, . . . , N. We obtain the following 



r (p+fc+") i 



,6 c=n T—n \ ' \ / 



r=p+| s=pT=0 



(2r-25)! (n! (2g)!! fs - v\ vr 

' {2r-2s-by.{p + b + ny.{2s-2p)U\ T )^ ' {n-p-2r+6+2.-2r, 2.-6- 

r=p+l+| 5=p T=0 ^ ^ ^ • 

(2r-2s)! n! (2g)!! 

' (2r - 2s - 6 - 1)! (p + 6 + n + 1)! (2s - 2p)!! 

^ J" ^ ("-'■) ^(n-p-2r+6+2s-2r+l,2r-6-2s+2T-l)- ('^•'^) 

case 3: h is odd and n = p. We have 
case 4: 6 is odd and n = p + 1, . . . , N. In this case we get 



6+1 



^ (p+b+n) j 



(2r-2s)! n! (2s)!! fs-p 



(2r - 2s - 6)! (p + 6 + n)! (2s - 2p)!! V T 

[ (p+b+n + l) ] 

"^(n-p-2r-+&+2s-2T, 2r--6-2s+2r) + ^ 2rH^ 

r=p+^ s=p T=0 

/p + b + n + l\ (_^y.+n(r\ (2r-2s)! ri! (2s)!! 

V 2r / (2r-2s-6-l)!(p + 6 + n+l)!(2s-2p)! 

( T ) ^"''"•^ "'^(n-p-2r+6+2s-2r+l,2r-6-2s+2T-l) (^■^) 
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Let us begin by considering the system f l4.2p for 6 = 0, so that it has an 
equal number of equations and unknowns. Moreover, the matrix Gm, n is non 
singular (because in the case Lp = and by the condition 1. of hyperbolicity, 
the system must have only the zero solution); as a result it yields 

YJX-^^ = 0, for n = p,...,N. (4.7) 

Now we have as an equation fl4.2p with b=l,...,N — pas well as (14. 7p . In 
particular, fl4.7p with n = p, for fl4.3p gives 

¥^,^^Xjoy>=0. (4.8) 

From this result we see that ^p^U^ = is one of the eigenvalues, because in 
this case one of the equations is an identity, so that we have less equations 
than unknowns to determine the eignvectors. If we look for other eigenvalues, 
that is ip^W^ 7^ 0, then (14. 8 p gives X'^qq-I^ = 0. Using this equation, together 
with dOD and we get 

Y,^^-'^^ = . (4.9) 

so that the term with n = p in the system (14. 2 p can be omitted. 

Summarizing the results obtained until now, we have found the eigenvalue 
(p^U^ = (which does not depend on Gm,n) and, for the other eigenvalues, 
the system 

^(0,0) 

Y,^X'^''=0, ioTn = p,...,N, 

En=p+iG'„.,„yo!n'^' = 0, for b=l,...N-p and m = p + b,...,N. 

(4.10) 



Let us now repeat the above steps, but with b = 1. In this case the third 
equation of f l4.1Up has again an equal number of equations and unknowns 

^71- 
b,p 

y;i^-> = 0, for n = p + l,...,N, (4.11) 



(because we have dropped Yfjl^""'^) and gives the solution 



and these equations replace the third equation of fl4.10p for b = 1. Now we 
note that by using (14. 4p and (14.60 . the second equations of (14.100 and (14. lip 
for n = p + 1 become 

^.U^j2^iP + 1)!(2P)!! ^(to7^ + ^j^iP +l)iP+ l)!(2j>)!!2 Xj-> = 0, 
-^12^(P + 1)2(P + mm ^(to^ - V^U'^(2^TW(P + + 1)!(2P)!! X^;^^ 
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This is a homogeneous system of two equations and two unknowns X^-^"^'^'' 
and X'l'^"^^'' , so that we have the following possibilities: 

• If the coefficient matrix is singular, we obtain the eigenvalues from the 
following: 

which allows us to verify that its eigenvalues do not depend on Gm,n, 

• Its solution is XJq)'^'' = and X^^^''-^^" = 0. 
From this second possibility it follows that 



n:p;? = o, (4.12) 

as seen from fl4.4p and (14 .6^ (Note that Y^Jl^""''' is a linear combination of 
■^{h k)^ with h + k = n — p; in our case h + k = 1; in other words, Yh^p^i^ 
is a linear combination of X^^-^"^^^ and X^^'^^'' which are zero in the present 
case). From f l4.12p it follows that in the third equation of fl4.10p we can 
omit the term with n = p + 1. Summarizing the results of this new step, 
we have found the set of eigenvalues (the solutions of ^p^U^ = and of 
{ip^U^y — 2^V^ = 0) and, to determine the other eigenvalues, the system 

xjH:ki'' = ^Mh+k<i 

Y,^X'^''=0, foTn = p+l,...,N, 

En=p+2Gm,nY,''n'''' = 0, for 6 = 2, . . . iV - p and m = p + 6, . . . , iV. 

(4.13) 

We also note that the set S"^ does not depend on Gm,n- Let us now iterate 
this procedure i] times and find 

• A set of eigenvalues not depending on Gm,n, 

• The following system, for the determination of other eventual eigenval- 
ues 

Xj-^^=0,for/i + A;<r7-l 

YJX'^^ = 0, for q = 0, ... ,1] — 1 and n = p + q, . . . , N , 
Etp+rj^m^nY^'X''" =0, foT b = ^ , . . . N - p ^nd ni = p + b, . . . , N . 

(4.14) 
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It remains to prove that, starting from this hypothesis, it follows that it holds 
also with ij+l instead of r]. The system given by the third equation of f l4.14p 
with b = r] has the solutions 

"^vX"^" = 0, forn = j9 + ?7, . . . , iV, (4.15) 

and these equations replace the third equation of (14.141) for b = rj. We also 
note that by using (14. 4 p and (14. 6p . the second equation of (I4.14p and (I4.15P 
for n = p + rj constitute a system of homogeneous equations in the unknowns 
■^{h k)^ with h + k = 1]. Thus we have an equal number (77 + 1) of equations 
and of unknowns and this system does not depend on Gm.,n- By imposing that 
its coefficient matrix is singular, we obtain some eigenvalues which, together 
with S"^, constitute the new set S^^^. If we look for other eigenvalues, then 
this system has only X'^^"^^ = with h + k = r] as a solution. On the other 

hand, this solution implies = 0, as seen from (14. 4p and (14. 6p . As a 

result, we can now omit the term with n = p + rj in the third equation of 
(I4.14p . In other words, we have found the set S^^^ and the system (I4.14p 
with T] + 1 instead of i]. This completes the proof of this property. 

We note that our system (14. 2 p has been gradually replaced by (14. 7p with 
n = p ai the first step, by the second equation of (I4.10p and (14. lip with 
n = p + 1 at the second step, and so on. In other words, our system (14. 2 p 
can be replaced by 

y,]p+v=0, for g = 0,...,r/, (4.16) 

and on this system we have to impose first the subsystem with 77 = to 
determine X^^-^'q^^ , then the subsystem with i] = Ito determine X^'^^'^J*' where 
h + k = 1, and so on for increasing values of r]. In particular, (I4.16p . for fixed 
T], will be a subsystem to determine Xj^"j^^ where h + k = i]. 

Thus we have proved that the wave speeds do not depend on Gm,n- This 
fact allows us to use (14. lap , without loss of generality. A natural question 
is the following: Will this choice satisfy the above condition on the non 
singularity of the matrix Gm,n for m, n = p + i], . . . N7 We prove that this 
is indeed the case. First, we see that (14. lap with the change of integration 
variable p = ^x, becomes 

T POO 1 /i\m+n+3 /-oo 

(4.17) 
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Integrating by parts the expression e ^x^ dx , we get 



/■oo poo poo 

/ e'^'x^ dx = \-e~'^x'^\^ + p / e'"" x^^'^ dx = p I e'^'x^'^ dx. 
Jo Jo Jo 



Iterating the integration other p — 1 times we arrive at 

POO 

e~^x^dx = p\ / dx = p\. 



As a result, (I4.17P becomes 



m+n+3 



[m + n + 2)\. 



It follows that 



Crp+r)+ l,p+ri Gp-^^r|+ l,p+r)+l 



G 



N,p+r) 



G 



N,p+ri+l 



Gp+n, N 

Gp+.q+l, N 



G 



N,N 



(2p+2r;+2)! ( k \ 



^ / , \ 2p+2»7+3 



fc2 



_ A 

e fc 



(2p+2r;+3)! / k 
k^ I 7 



jp+ri+N +2)1 ( k\ 



2p+2»7+4 ^ 



p+ri+N+S ^ 



{p+ri+N+2y. 

(p+rj+N+3y. 
fc2 
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_A\ N-p-r^+l 

e k 



N+l—p—r) 

{2p + ri + i) 

1=3 



{2p + 2ri + 2)l (2p + 2ry + 3)!^ • 
{2p + 2r] + 3)\ (2p + 2?7 + 4)!^ • 

(p + r/ + 7V + 2)! {p + 2ri + N + 3)\^ ■ 

J2 (2p + r? + 



e k 



1=3 



(2p + 2?7 + 2)! (2p + 2r; + 3)! 
(2p + 2r/ + 3)! (2p + 2r/ + 4)! 

{p + rj + N + 2)1 {p + 2r] + N + 3)\ 

Let us now consider the determinant 

a! (a + 1)! 

(a+1)! (a + 2)! 



• (p + 7/ + iV + 2)!(^ 

• (p + ^ + iV + 3)!(^) 

(2iV + 2)!(^) 

N—p—rj 

+ E ^ 



{p + r] + N + 2)\ 
{p + r] + N + 3)\ 

{2N + 2)\ 



N—p—r) 



N—p—rj 



N—p—r] 



(4.18a) 



D 



a, d 



(a + c)! (a + c+1)! 
{a + d)\ {a + d+iy. 



{a + b)\ 
(a + fo + 1)! 

(a + 6 + c)! 

{a + b + d)\ 



{a + d)\ 
{a + d+iy. 

{a + d + c)\ 

{a + 2d)\ 



and let us sum to the hne beginning with (a + c)! the previous one multiphed 
by —(a + c). We obtain 



Da. 



a! (a + 1)! 

(a+1)! 

(a + c)! 

{a + d)\ 



{a + b)\ 
b{a + b)\ 

b{a + b + c-l)\ 

b{a + b + d-l)\ 



{a + d)\ 
d(a + d)\ 

d{a + d + c- 1)! 

d{a + 2d~ 1)! 
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where we have taken into account that {a + b + c)\ — {a + c){a + b + c — 1)\ 



[a + b + c - l)\{a + b + c - a - c) = b{a + b + c - 1) 



Da,d = a\ d\ 



[a + l)\ 
{a + c)\ 



{a + d)\ 
a! d\ Da+i,d^i- 



{a + b)\ 
(a + 6 + c- 1)! 
{a + b + d-l)\ 



It follows that 
{a + d)\ 

[a + d + c- 1)! 

{a + 2d- 1)! 



By iterating the procedure r times, we obtain 

Da,d = a\{a+l)\---{a-l + r)! d\ {d - 1)\ ■ ■ ■ {d + 1 - r)\ Da+r,d-r. 

As a result, for r = d this last expression becomes 

Da,d = a\{a + l)\---{a-l + d)\ d\ (d - 1)! ■ ■ • 2! |(a + d)\\ 
= a\ {a + l)\---{a + d)\d\ (d- l)!---2! 

By applying this result with a = 2p + 2rj + 2, d = N —p — t], we see that the 
determinant in fl4.18p is equal to 

{2p+2r]+2)\ (2p+2?7+3)! ■ ■ ■ {p+r]+N+2)\ {N-p-r])\ {N-p-r]-l)\ ■ ■ ■ 2! > . 
It also follows that the matrix Gm,n for m, n = 0, . . . , N is positive definite. 
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